Accounting for Breit interaction in actinide and superheavy element compounds: 

1. General remarks. 
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The incorporation of the Breit interaction in atomic and molecular calculations is discussed in 
the framework of four-component all-electron and two-component relativistic effective core potential 
(RECP) formalisms. Contributions of the Breit interaction between different core and valence shells 
are studied in the Dirac-Fock approximation for the uranium, plutonium, ekathallium (^=113), 
ekalead (Z=114) and other heavy atoms. It is shown that the two-electron Breit effects between the 
valence electrons can be neglected for "chemical accuracy" (I kcal/mol) of calculation of spectro- 
scopic properties of systems containing superheavy elements and actinides whereas large core-core 
and core-valence Breit contributions can be efficiently described by one-electron RECP operators. 
Different versions of the generalized RECPs with the Breit interaction taken into account are con- 
structed for uranium and plutonium and tested in comparison with the corresponding all-electron 
four-component calculations. 



PACS numbers: 31.10.-l-z, 31.15.Fx, 31.15.Ne, 31.30.Jv 



Introduction. 

The Dirac Hamiltonian which accounts for relativistic effects for an electron in an external potential v can be 
written in the form 

hP ~ cd-p + (3mc^ + V , (1) 

where a and P are Dirac matrices, c is velocity of light. The generalization of this equation to the many-electron 
case is not as straightforward as is in the nonrelativistic case. For few-electron systems, the quantum electrodynamics 
(QED) calculations can be usually performed However, the perturbative QED approach is not applicable 

in practice to many-electron systems such as heavy atoms and their compounds when correlation effects are also 
very important. Instead, the Hamiltonian formalism with some approximations for one- and two-electron parts is 
employed. For many-electron systems, the effective Hamiltonian (in which relatively small irreducible three-electron 
etc. interactions are neglected) is usually written as 

H = E'^? + EG^^' (2) 

i i>j 

where Gij is an effective operator which describes interaction between electrons i and j. It can be most naturally 
derived from QED. In the first QED order, this operator describes interelectronic exchange by a single virtual photon. 
The problem is that the Gij operator is gauge dependent in this approximation, e.g., in Feynman gauge 

Gfj i^ij ) - (1 - a» -aj ) '-^ , (3) 

rij 

where is the frequency of the exchanged photon divided by the velocity of light, and in Coulomb gauge 



di-Vi, 



^ - cos{uJij rij)-l 

'^j'^j^ , ,2 



(4) 
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As is shown in |3'|, operators and Q give equal mean values with one-configuration wave function constructed 
of one-particle solutions of eq. ^ if v{'r) is a local potential. However, it is demonstrated both theoretically and 
computationally in 0, 0] that the result is gauge-dependent in a general case. Then, it was clarified in0 that the 
gauge dependence remains if one goes beyond the no-virtual-pair approximation as well. It was noted in 7] that when 
solving the eigenvalue problem with Hamiltonian (0) , only so called ladder diagrams of a multi-photon exchange are 
taken into account. When two-photon crossed diagram is also considered, the two gauges give the same result in the 
leading relativistic order, 0{a^) , where q;w1/137 is the fine structure constant. Because the crossed diagram does 
not contribute in this order in the Coulomb gauge 0, the electron-electron interaction in form Q is usually more 
appropriate for electronic structure calculations based on eq. 101. 

In this paper we study the Breit interaction in the Dirac-Fock approximation. For the reason discussed below, we 
will further use the low-frequency limit of eq. I0J 

GfjiO) = l/rij - di-dj/nj +^ [di■dJ-{d^■nj){dynJ) /r^,j\ jr^ . (5) 

The first term in eq. (O describes the instantaneous Coulomb interaction between electrons and the following two 
terms correspond to the original Breit interaction (BI) . The second term in eq. describes instantaneous magnetic or 
Gaunt interaction which is usually a dominant part of BI for an atomic system and the third term describes classical 
retardation of the Coulomb interaction between electrons. Operator ||2Jl with electron-electron interaction (O is called 
Dirac-Coulomb-Breit (DCB) Hamiltonian. If the retardation part or both the retardation and magnetic parts of BI 
are neglected, one goes to the Dirac-Coulomb-Gaunt (DCG) or Dirac-Coulomb (DC) Hamiltonians, correspondingly. 

The terms of order of O{ojfjrij) , which represent higher order retardation, are omitted in low- frequency approxi- 
mation ||SJ) for electron-electron interaction I^J . This approximation is not valid for BI between negative (positronic) 
and positive (electronic) energy states, when cuijc w 2m(? . When low-frequency approximation (jsj is applied to derive 
the positive and negative energy one-particle states as solutions of the Dirac-Fock-Breit equations, the errors for the 
energies are of order of 0{a^) as is shown in ,8]. A serious disadvantage of the high-order retardation terms is that 
they lead to the non-Hermitian exchange part of the Dirac-Fock potential 0. It is pointed out in that inclusion 
of higher order retardation effects in calculations without self-energy terms is not reasonable because, at least, part 
of these effects cancels each other. From the computational point of view, the higher order retardation terms lead to 
only a small correction to the cumulative BI effect (e.g., see, ^^i^, if necessary, they can be taken into account 
posteriori as a first-order perturbation. 



I. RESULTS OF ATOMIC DC(B)-CALCULATIONS. 

In Tables lll < VII we report the results of calculations for a series of heavy elements by the self-consistent field (SCF) 
method with the DC, DCG and DCB Hamiltonians. The finite-difference (numerical) HFD code llj was developed 
to take into account the Gaunt and Breit interactions. The absolute values of one-electron energies for most of 
spinors of the uranium atom are decreased approximately on 0.2-0.5% when BI is taken into account in the SCF 
procedure. The largest change in one-electron energies among the spinors with the same principal quantum numbers 
is observed for the j)i/2 states (see Table iHlj) . The BI contribution to the one-electron energies has the opposite 
sign for the 6(i, 4/, 5/ spinors as compared to those for other spinors. The energies of the 5/5/2 and 5/7/2 states 
whose occupation numbers are changed in some of transitions between the low-lying states of the uranium atom 
and its compounds are decreased by more than 1%, on 820 and 1000 cm"-'^, correspondingly. Similar behavior was 
observed in Q for the 4/ and bd spinors of the mercury atom. The inversion of the sign is caused by the change 
in the electric part (Coulomb plus exchange) of the Dirac-Fock potential when spinors are relaxed due to the Breit 
interaction |l[l3. The differ ence between matrix elements of the hyperfine interaction in calculations with the DC, 
DCG and DCB Hamiltonians takes place only due to the spinor relaxation. The retardation corrections give about 
10% of the total Breit contributions to the one-electron energies. They can achieve more than 30% for the hyperfine 
interaction and usually have the opposite sign as compared to the Gaunt corrections. 

Table m demonstrates that the uranium atom has low-lying states with different occupation numbers of the 5/ shell. 
Accounting for the BI contribution is especially important for transitions with the excitation of 5/ electrons because 
of the great contribution of BI to its one-electron energies. The BI contribution to the transition energies is up to 
3000 cm~^(and sometimes is comparable with the DC transition energy). Similar situation is observed for plutonium 
(Table 1^) and is expected for other actinides. Obviously, the cumulative effect caused by BI for other calculated 
atoms grows for elements of the same group with increasing the nuclear charge Z and can achieve few hundreds wave 
numbers for very heavy elements (see Tables^^^nj- Nevertheless, the BI between valence electrons is decreased with 
growth of Z mainly because of enlarging the average radii of the valence shells as is shown below. 
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II. DISADVANTAGES OF DC(B)-BASED CALCULATIONS 

In the previous section, atomic finite-difFerence SCF calculations based on the DC, DCG and DCB Hamiltonians 
are discussed. From the formal point of view, correlation calculations with these Hamiltonians can provide a very 
high accuracy of physical and chemical properties for molecules containing heavy atoms. However, such calculations 
are not widely used yet for such systems because of the following theoretical and technical complications 

1 Too many electrons are treated explicitly in heavy-atom systems and too large basis set of Gaussians is required 
for accurate description of the large number of oscillations, which valence spinors have in the presence of heavy 
atoms. Because the number of two-electron integrals grows as the fourth degree of the number of basis functions, 
it leads to serious complications of their calculations and, especially, of their transformation to a molecular basis. 

2 When the Dirac formalism with four-component spinors is used, the number of two-electron integrals to be 
computed is strongly increased as compared to the one-component (nonrclativistic or scalar-relativistic) and 
two-component (spin-dependent) cases: 

(a) The number of kinetically-balanced two-component ( "2c" ) uncontracted basis spinors required for descrip- 
tion of the S'mall components, Ng'^, can be estimated as 2iV^'^, where N'l'' is the number of basis spinors 
for the Large components; so the total number of the relativistic four-component ("4c") basis spinors 
^bas ^ ■^-^L^ ^^"i t^*^ number of calculated two-electron integrals is 

Ntint - (l+2-22+24)iV|,-,„, ^ 25.iV|,=,„, , 

where N^^j^^ is the number of two-electron integrals with only the large components when the Coulomb 
electron-electron interaction is taken into account. Even more computational effort is required when the 
magnetic and retardation electron-electron interactions are taken into account. To reduce the number of 
two-electron integrals, different two-component relativistic approximations are actively developed first of 
all to be used in molecular calculations (e.g., see 0). 

(b) The number of basis 2c-spinors, iV^^g, is twice more than the number of nonrclativistic basis one-component 
("Ic") orbitals, Nl^^, therefore 

(the division by factor of two appears due to different symmetry properties of two-electron integrals for Ic 
and 2c (or 4c) cases) and A^|e/„t ~ 4-;-6Af|g/^, where A^|e/,5? is the number of two-electron integrals in the 
spin-orbit basis set which are required to be saved in the computer memory. 

3 Eq. Q as well as the Dirac-Fock equations has both negative- (positronic) and positive-energy (electronic) solu- 
tions. This circumstance leads to so called "continuum dissolution" when eq. (0) is applied to electronic structure 
calculations 16j. To prevent it, one can use only the positive-energy solutions for constructing many-electron 
basis functions ("no-virtual-pair" approximation). The way to go beyond the no- virtual-pair approximation 
is to use the normal ordered second-quantized representation of eq. Q with respect to the "vacuum" level in 
which all the negative-energy states are occupied. Evaluation of the Dirac-Fock(-Breit) bispinors is required 
for separating the negative-energy and positive-energy solutions in the one-particle basis. Because the virtual 
Dirac-Fock(-Breit) bispinors are not optimal for precise correlation calculations of molecules, some additional 
efforts can be required to generate an appropriate set of spinors for such calculations. 

To avoid complications described in items 2a and 3, small components can be excluded explicitly from calculations 
llH. To reduce also the difficulties described in items 1 and 26, Hamiltonian can be replaced by an effective 
Hamiltonian 



in practical applications when inner core electrons do not play an active role and, besides, the original valence spinors 
can be smoothed in the heavy-atom cores. Hamiltonian © is written only for valence or "valence-extended" (VE) 
subspace of electrons denoted by indices i„ and jy (in the VE case, some outermost core shells are also treated 
explicitly). U^* is an relativistic effective core potential (RECP) or relativistic pseudopotential (RPP) operator that 
is usually written in the radially-local (semi- local) or separable form (e.g., see and references). It simulates. 
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in particular, interactions (which can include BI, see below) of the explicitly treated VE electrons with those which 
are excluded from the RECP calculations. Besides, the generalized RECP (GRECP) operator 0,113 can be used, 
which includes both the radially-local and separable terms (additionally, the GRECP method involves possibilities 
of the Huzinaga-typc pseudopotential and other advantages, see below). In Eq. h^'^^'^ is the one-electron 
Schrodinger Hamiltonian 

jjSchr ^ _ ly2 _^ ^nuc ^ j^y-j 

Contrary to the four-component wave function used in DC(B) calculations, the pseudo-wave function in the (G)RECP 
case can be both two- and one-component. The use of Hamiltonian (|SJ| instead of |(2J) is, of course, an additional ap- 
proximation and the question about its accuracy appears. As was shown both theoretically and in many calculations, 
the usual accuracy of the radially- local RECP versions is within 1000-3000 cm~^for transition energies and an ap- 
propriate level of accuracy is expected for other chemical and physical properties. 

III. GENERALIZED RELATIVISTIC EFFECTIVE CORE POTENTIAL. 

In a series of papers (see 0, |23| and references) , we have introduced and developed the GRECP concept which 
allows one to attain practically any desired accuracy, while requiring moderate computational efforts. 
The main steps of the scheme of generation of the GRECP version with the separable correction are: 

1. The numerical all-electron relativistic calculation of a generator state is carried out for an atom under consid- 
eration. For this purpose, we use the atomic HFDB code [TllIT^. 

2. The numerical pseudospinors (pnij{r) are constructed of the large components ipniji''') of the outer-core (C) and 
valence (V) HFDB spinors so that the innermost pseudospinors of them (for each I and j) are nodeless, the next 
pseudospinors have one node, and so forth. These pseudospinors satisfy the following conditions: 

Z = 0,1,...,L,^ j = \l±^\, 

71 — 71^, 71^, . . . , Tly^ 

where L is one more than the highest orbital angular momentum of the inner core (IC) spinors. The leading 
power 7 in the polynomial is typically chosen to be close to L -I- 1 in order to ensure a sufficient ejection of 
the valence and outer-core electrons from the IC region. The Ui coefficients are determined by the following 
requirements: 

• {^nij} set is orthonormalized. 



y and its first four derivatives match i/Jnij and its derivatives, 
y is a smooth and nodeless function, and 

ipnij ensures a sufficiently smooth shape of the corresponding potential. 



Rc is chosen near the extremum of the spinor so that the corresponding pseudospinor has the defined above 
number of nodes. In practice, the Rc radii for the different spinors should be chosen close to each other to 
generate smooth potentials. 

3. The Unij potentials are derived for each /=0, . . . , L and ± ^| for the valence and outer-core pseudospinors 
so that the ipnij are solutions of the nonrelativistic-type Hartree-Fock equations in the jj-coupling scheme for a 
"pseudoatom" with the removed IC electrons. 




(9) 
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where Zic is the charge of the nucleus decreased by the number of IC electrons, J and K are Coulomb and 
exchange operators calculated with the ^nij pseudospinors, Snij are the one-electron energies of the corresponding 
spinors, and en'nij are off-diagonal Lagrange multipliers (which are, in general, slightly different for the original 
bispinors and pseudospinors). 

4. The GRECP operator with the separable correction written in the spinor representation is as 
U = f/„„Lj(r) + E E { [UnjAr) - Un^Lj{r)]Pij 

l=0j=\l-l/2\ 

+ E [UnMir) ~ Unjjir)] P„^i, -|- E Pnjj [Unjj{r) - Unj,{r)] (10) 

ric Tie 

~ E/ -^"c'j 

where 

i _ j 

Py ^ E \ljm){ljm\, Pnjj = E \njj^^){njjm\, 

rn——j m—~j 

\ljm'^(^ljm\ is the projector on the two-component spin-angular function xijm, \ncljm^(^ncljm\ is the projector 
on the outer core pseudospinors ^njjXijm, and J = L + 1/2. 

Two of the major features of the GRECP version with the separable correction described here are generating of the 
effective potential components for pseudospinors which may have nodes, and addition of non-local separable terms 
with projectors on the outer core pseudospinors (the second and third lines in eq. <^lQ\i ) to the standard semi-local 
RECP operator (the first line in eq. (|1U|) ). Some other GRECP versions are described and discussed in papers 
[l9, 22] in details. 

The GRECP operator in spinor representation (|10|l is mainly used in our atomic calculations. The spin-orbit 
representation of this operator which can be found in 19] is more efficient in practice being applied to molecular 
calculations. Despite the complexity of expression (|1L)|) for the GRECP operator, the calculation of its one-electron 
integrals is not significantly more expensive than that for the case of the standard radially-local RECP operator. 



n'Jj 



IV. ACCOUNTING FOR BREIT EFFECTS IN GRECP 



Let us consider the contributions of BI between electrons from different shells to the total energy of a heavy 
atom 122^. The following estimate can be applied (e.g., see p^ ]! 

{a,b\{di-dj)\a,b) ^ ^{{va-Vb)) . 



For an uncoupled one-electron state a one has 

{a\d\a) 



ma\ 



where Z* is an effective charge of the core that is experienced by the electron in the a-th state, Z*=Z— iV", Z is 
the nuclear charge, TV" is the number of core electrons with respect to the a-th state (the core radius i?" ^ {a\r\a)). 
Besides, (l/r'12) can be estimated as (l/r) for the outermost (from a and h) one-electron state: 



1 



{ab\ \ab) ^ min 

ri2 



(a|i|a),(fo|i|6) 



min[Z*,Zb] . 

Thus, the absolute value of BI between electrons in states a and b can be roughly estimated as 



\Bab\-^a'z:z:- mm [Z:,Zt] 



(11) 
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When eq. @ is applied instead of eq. (0), it means neglecting (or, in some sense, freezing) the BI between the VE 
electrons. In the case of a=v and b=v', one has Z* and 

\Byy, \ ^a^ K. ^•10^'' a.u. w 10 cm"^ . 

Thus, this contribution is negligible for the "chemical accuracy" (about 350 cm^^) of calculations. This is demon- 
strated in Tables IVIII and IVIIII where we calculated the BI contributions to the total and transition energies of the 
uranium atom. In these calculations we have subdivided all the electrons on the "core" and "valence" subspaces by 
different ways and the BI contributions between the valence electrons are neglected for each of such subdivisions (i.e., 
only core-core and core-valence contributions are taken into account). We have also calculated the BI contributions to 
the total energy between the valence- valence, core-valence and core-core shells and have estimated them for different 
s shells applying eq. The results are collected in Table Hxl In this table we have calculated Z* as {a\\/r\a) and 

have introduced the normalizing factors Cn = 0.34 as some average multiplier for the BI contributions ("direct" plus 
"exchange") per an electron pair within the considered s shells (see the footnote to Table Hxl for more details). It 
may be seen from Tabic HXl that the estimates from agree with the calculated values mainly within the factor of 
two. The suggested estimate is expected to be substantially cruder for higher total electronic momenta. 

The valence-valence contribution is usually decreased with the increase of Z . There are several reasons for such 
behaviour. The expression for BI contains double integration on radial variables of the products of large and small 
radial components of Dirac bispinor 



j dr, I dr2^i(ri)V'c(''2)^^fe(ri,r2)V'6^(ri)V'd(r2)r?r2' , (12) 



where Uk{ri,r2) — f^/r^^, r'< = min[ri,r2] and r> — max[ri,r2]. For high Z radial parts of the large and small 
components of the valence spinors have a number of oscillations, which usually do not coincide. It leads to oscillations 
in the products '0^(r)-'!/;^(r)-r^ (see Fig.^ and large cancellation in integral H12|l as compared to the case of integrating 
over the amplitude of the integrant (see Table llX)l : 



j dr, j dr2 I i>i{r,)^'Ar2)uu{r,,r2)rb{ri)rd{r2) \ rlrl ■ (13) 



In particular, it also means that the "incompensated" BI contribution between the valence electrons is mainly taken in 
the valence region. As one can see from Fig. ^ the absolute value of the ril}\ (r) • ripl (i^) is decreased with the increase 
of Z. It may be easily understood when the kinetic balance condition is applied to derive the small components from 
the large ones in the valence region: 



r-r{r) = f + [r-i^'ir)] . (14) 

The atomic radius is increased with the increase of Z and, thus, the large components of the valence spinors become 
smaller and smoother. It leads to the decrease of the small components in accord to eq. (|14|l . 

The core electrons can be considered as frozen when studying majority of physical-chemical properties and processes 
of practical interest. Since they usually form closed shells in calculations, the total contribution of the direct part of 
BI between the core and valence electrons is equal to zero and only the exchange part of this interaction gives nonzero 
contribution. For the innermost Is-shell (a=c) of a heavy atom, one has Z*^Z'^100 by the order of magnitude and 

\Bcv\ - a^-100 « ^-10"^ a.u. « 1000 cm^^ , 

that is quite essential for calculations on the level of "chemical accuracy". (We found in the present calculations 
that contribution from the exchange interaction is usually larger than the direct one in the cases of pairs of spinors 
{nsi/2,msi/2); {nsi/2,mpi/2); (n-Si/2, mp3/2); inpi/2,mpi/2); (npi/2, 771^3/2) and smaller in the cases of pairs of 
spinors (nsi/2, md3/2); (nsi/2, mds/s); (nsi/2, TO/5/2); (nsi/2, TO/7/2); (npi/2, TOP3/2); {npi/2,ind^/2)\ ('^Pi/2, '^1/7/2); 
(np3/2, TO(i3/2); (n-P3/2j TO/5/2); («P3/2, TO/7/2).) Therefore, one can just take into account the Bcv contributions (and 
neglect the B^yi terms) when generating (G)RECPs with only valence electrons treated explicitly, where index c 
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runs over the core electrons only. The latter are explicitly excluded from the (G)RECP calculations, therefore, the 
effective operator for B^v acting on the valence shells has the same spin-angular dependence as the conventional 
radially-local RECP has. (For higher accuracy, GRECP can be used, which, besides, can efficiently treat other 
QED effects even for highly charged ions of heavy-element systems.) Thus, BI can be taken into account directly 
when the HFDB calculations are performed (see item 1 in the previous section) to generate valence bispinors and 
their one-electron energies, but in the procedure of inversion of the two-component HE equations for generating the 
(G)RECP components (see item 3 in the previous section), the conventional interelectronic Coulomb interaction can 
be used instead of the Coulomb-Breit one. Afterwards, one should consider only the Coulomb interaction between 
the explicitly treated electrons in the correlation (G)RECP calculations. In principle, the errors due to neglecting BI 
between the VE electrons (as well as the errors in reproducing the interelectronic Coulomb interaction due to spinor 
smoothing) can be compensated with good accuracy by the term-splitting GRECP correction 19, 20] . Alternatively, 
the procedure of restoration of proper shapes for the four-component VE spinors in the heavy-atom cores p3 , l25l | 
after the molecular GRECP calculation can be applied, then both BI between the VE electrons and the proper 
interelectronic Coulomb interaction are directly taken into account with the restored wave function. 

The Breit interaction for both the core electrons. Bed, can be of the same order of magnitude as the Coulomb 
interaction between them. However, Bed does not contribute to "differential" (valence) properties directly, whereas 
its indirect contribution can be easily taken into account with the help of the (G)RECP technique. 



V. RESULTS OF THE GRECP CALCULATIONS 



In the given section we consider accounting for BI by the GRECP method only for actinides for which the Breit 
effects can attain a few thousands wave numbers even to the energies of lowest-lying electronic excitations. Accounting 
for BI in different RECP and GRECP versions for superheavy elements and actinides will be considered in details 
in our forthcoming paper . For the uranium and plutonium atoms we constructed different 24- and 26-electrons 
GRECPs, correspondingly, which effectively account for the Breit effects. The GRECP(/3) and GRECP(/2) for the 
uranium atom were constructed using the generator states with the occupation numbers of the 5/-shell equal to 3 
and 2. In the case of the plutonium atom the configurations with the occupation numbers of the 5/-shell equal to 6 
and 5 were used, respectively. The results of the GRECP/SCF calculations as compared to those of all-electron HFDB 
calculations are presented in Tables HI and UTI 

The errors of calculations with these GRECPs can be collected into two groups. The errors for transitions without 
the change in the occupation number of the 5/ shell are rather small. The errors for transitions with the change in 
the occupation number of this shell can achieve few thousand wave numbers. The latter errors have a systematic 
nature and are mainly connected with two facts: (a) the 5/ shell of uranium and plutonium is described with the help 
of nodeless pseudospinors in the present GRECP versions and (b) the relaxation of the [Xe4/^^]-like core shells is 
not taken into account since they are excluded from the GRECP calculations explicitly. These errors can be reduced 
significantly if one includes the 4/^'*5s^5p^ electrons in the calculations for uranium and plutonium explicitly. Another 
way is to apply the self-consistent (SfC) GRECP version [l9ll27j. In the SfC GRECP case, the circumstance is taken 
into account that the GRECPs generated for different occupation numbers of the outermost core shell(s) {N^f'^ in 
our case) are somewhat distinguished due to above mentioned reasons (see 19, 22] for more details). The dependence 
on N^j'^ is introduced explicitly to the SfC GRECP operators, which are constructed for uranium and plutonium on 
the basis of the above generated GRECPs. As it can be seen from Tables III and ITU the SfC GRECPs allow one to 
increase the accuracy of the calculations without extension of the space of explicitly treated electrons up to an order 
of magnitude and more. 
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TABLE I: Energies of transitions between the states averaged over nonrelativistic configurations of the uranium atom calculated 
by the DHF method with the Coulomb, Coulomb-Gaunt and Coulomb-Breit two-electrons interactions taken into account and 
by the SCF method with diflferent GRECPs (in cm-^). 
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36227 
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57 


-4 


5f7s^6d^7p^ 


17214 


-11 


-3 


-29 


-81 


-29 
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42373 


-47 


2 


-15 


-29 


-15 




54714 


-147 


9 


-11 


54 


-9 


5f7s''6d^ 
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-736 
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-424 
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15053 


-44 


1 


34 


63 


-3 
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14953 


-19 


-2 


-68 


-112 


-25 




38863 


-54 


4 


-33 


-43 


-21 


5f6d^ 


33874 


-84 


2 


63 


109 


-1 


5f6d'7p' 


32194 


-80 


1 


3 


11 


-11 
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-114 


6 


26 


65 


-15 
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-20 
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-17 
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-13 


-5 
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42710 


91 


-8 
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23 


32 
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-113 


9 


-103 


-20 


-24 




19232 


-11 


-4 


28 


-28 


-27 


5f7s'6d^ 


45402 


-52 


2 


-4 


-17 


-15 


5f&d^ 


54780 


-180 


11 


-115 


-27 


-32 


5f7s^6d^ 














5f7s^6d^ 


29773 


1964 


-290 


2728 


1285 


-141 


5f7s''6d^ 














5f7s^6d^7p^ 


18189 


156 


-19 


432 


242 


82 


5f7s^6d^ 


49233 


106 


-10 


340 


204 


99 




7455 


-134 


11 


-256 


-157 


-80 


5f7s^6d^7p^ 


21056 


-13 


-5 


93 


37 


-17 


5f7s^6d^ 


48058 


-59 


2 


-7 


-15 


-24 


5f6d'^ 


54003 


-210 


14 


-284 


-181 


-100 


5f7s^6d^ 














5f 


100587 


-1313 


187 


-534 


392 


14 
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TABLE II: Energies of transitions between the states averaged over nonrelativistic configurations of the plutonium atom 
calculated by the DHF method with the Coulomb, Coulomb-Gaunt and Coulomb-Breit two-electrons interactions taken into 
account and by the SCF method with different GRECPs (in cm~^). 







Gaunt 


Retard. 


Absolute 


error of 26e-GRECPs 






contr. 


contr. 


(GRECP(/'^) (GRECP(/'^) (SfC GRECP 


Transition 


DC (DCG-DC) (DCB-DCG) 


-DCB) 


-DCB) 


-DCB) 


5f7s-' ^ 
















17217 


-56 


3 


12 


32 


10 




15697 


-16 


-3 


-27 


-57 


-26 


5f7s' 


39900 


-50 


3 


-16 


-27 


-15 




56908 


-120 


6 


8 


31 


8 


5f7p' 


66747 


-70 





-32 


-86 


-31 


5f7s^ ^ 
















44194 


-587 


84 


-450 


-119 


-70 


5f7s' - 














5f6d} 


19931 


-58 


4 





9 


-4 


5f7p' 


14878 


-68 


6 


-29 


-26 


-17 


5f 


35053 


-106 


10 


-23 





-24 


5f7s^ ^ 














5f7s'^6d^ 


-3803 


826 


-122 


960 


497 


6 


5f7s'^6d^ 














bf7s'7p' 


6650 


105 


-12 


72 


-26 


-25 


5f7s^6d^ 


15125 


-87 


6 


-42 


-3 


-4 


5f7s^6d^7p^ 


18262 


-13 


-3 


19 


-14 


-14 


5f7s^ 


35849 


65 


-4 


74 


3 


4 


5f7s^6d^ 


43812 


-51 


3 


3 


-7 


-6 


5f7s^ ^ 
















15880 


1812 


-267 


2585 


1580 


-124 


5f7s^6d^ 














bf7s^Gd^7p^ 


12319 


130 


-15 


273 


164 


63 


5f7s^6d^ 


12325 


-112 


8 


-146 


-94 


-49 


5f7s'6d^7p'' 


20424 


-14 


-5 


75 


42 


8 


5f7s^6d' 


42764 


84 


-7 


219 


140 


71 




47004 


-57 


2 


11 


3 


-5 


5f7s^ ^ 
















60152 


2927 


-431 


5022 


3412 


-411 


5f7s^6d^ 














5f7s^6d'' 


9050 


-134 


10 


-295 


-235 


-136 


5f7s^6d^7p' 


18114 


151 


-18 


533 


416 


208 


hf7s^&d^7p^ 


22344 


-16 


-5 


140 


107 


34 


5f6d^ 


24372 


-251 


19 


-522 


-413 


-230 


bf7s'Qd' 


49588 


97 


-8 


399 


314 


174 


5f7s^6d'^ 


49757 


-65 


2 


7 


1 


-15 
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TABLE III: Diagonal matrix elements of the hyperfine interaction (in MHz) and one-electron energies (in a.u.) obtained by 
the Dirac-Hartree-Fock method with the Coulomb, Coulomb-Gaunt and Coulomb-Breit two-electron interactions taken into 
account in calculation of the bf^ls^Gd} configuration of the ^^^U atom. 





Hyperfine interaction 


One-electron energy 




spinor 


DC 


DCG 


DCB 


DC 


DCG 


DCB 






-DC 


-DCG 




-DC 


-DCG 


Is 


-36798227 


133035 


7027 


-4279.2641972 


19.1907227 


-1.3823594 


2s 


-5773435 


14022 


-4873 


-806.1661987 


2.3414060 


-.2158518 


3s 


-1355941 


4565 


-1898 


-206.6114745 


.4637954 


-.0453239 


4s 


-364146 


1568 


-613 


-54.3244989 


.1054092 


-.0110359 


5s 


-89193.1 


406.9 


-151.3 


-12.5963766 


.0217699 


-.0021773 


6s 


-17527.1 


81.0 


-28.8 


-2.1379683 


.0038303 


-.0003502 


7s 


-1595.43 


5.72 


-2.16 


-.2023518 


.0002373 


-.0000114 


2iDi /o 


-1967173 


20229 


-1243 


-776.3775115 


3.9620074 


-.3915975 


3j3l/2 


-473306 


4570 


-516 


-193.1006247 


.7876421 


-.0776144 


4pi/2 


-124395 


1263 


-171 


-48.1943130 


.1833323 


-.0187803 


5Pl/2 


-28600.0 


296.6 


-39.4 


-10.1279114 


.0375639 


-.0037016 


6pi/2 


-4887.86 


55.87 


-7.01 


-1.3430981 


.0059178 


-.0005425 


2p3/2 


-136925 


1266 


-246 


-635.5705743 


2.5069539 


-.3943068 


3p3/2 


-33035.1 


275.9 


-69.4 


-160.3225666 


.4802565 


-.0777980 


4p3/2 


-8678.45 


74.70 


-20.67 


-39.5412400 


.0999325 


-.0177798 


5P3/2 




iO.o / 


-4. to 


Q nno'7/1 A 1 
— 8.U9z(44i 


.Ui / /9Z0 


-.0035444 


fiyjo /o 


-300.545 


2.773 


-.767 


-.9846770 


.0022487 


-.0005017 


34/2 


-17399.0 


94.4 


-11.9 


-139.0184571 


.3490051 


-.0371458 


44/2 


-4214.13 


22.32 


-3.95 


-29.7336816 


.0559536 


-.0058241 


54/2 


-780.353 


3.686 


-.691 


—4.3513548 


.0051781 


-.0004503 


64/2 


-49.0723 


.1622 


-.0351 


-.1927517 


-.0002264 


nnnn997 


34/2 


-6517.26 


29.25 


-5.78 


-132.4187434 


.2429070 


-.0357106 


44/2 


-1568.65 


6.29 


-1.75 


-28.1305532 


.0306448 


-.0055668 


54/2 


-286.521 


.910 


-.306 


-4.0416478 


.0004606 


-.0004474 


64/2 


-16.9990 


.0176 


-.0158 


-.1832747 


-.0003916 


.0000202 


4/5/2 


-773.708 


1.446 


-.190 


-15.2045521 


-.0051455 


.0026720 


5/5/2 


-62.1555 


-.1984 


.0271 


-.3470475 


-.0044479 


.0007116 


4/7/2 


-410.285 


.465 


-.143 


-14.7923447 


-.0166711 


.0026479 


5/7/2 


-32.0892 


-.1674 


.0122 


-.3195979 


-.0052414 


.0006863 



[26] N. S. Mosyagiri et al., in preparation. 

[27] A. V. Titov and N. S. Mosyagin, Structural Chem. 6, 317 (1995). 



11 



TABLE IV: Transition energies between states with ns^np^ configuration of the Tin {Z = 50), Lead {Z = 82) and Ekalead 
{Z — 114) atoms calculated by Dirac-Hartree-Fock method with Coulomb and Coulomb-Gaunt two-electron interactions. All 
values (excepting the relative differences) are in cm~^. 



Tin atom 


configuration, 


DC DCG 


absolute relative 


term 




difference difference (%) 



{5si/.^bi/y2),J = 3113 3153 40 L3 

(5si/25p;/25pJ/2)>>^ = 1 

(5s?/25p}/25P3/2).-^ = 2 5143 5139 -4 -0.1 

(5s?/25p^/2)> = 2 5941 5893 -48 -0.8 

(5s?/25p3/2), J = 15873 15820 -53 -0.3 



Lead atom 


configuration, 


DC DCG 


absolute relative 


term 




difference difference (%) 



{6si^2^pi/2),j = 

(6si/26pJ/26p3/2).-^ = 1 4752 4644 -108 -2.3 

(6s2/26pJ/26P;V2), J = 2 9625 9514 -111 -1.2 

(6S1/26P3/2): = 2 18826 18592 -234 -1.2 

(6s?/26P3/2), = 28239 27995 -244 -0.9 



Ekaload atom 


configuration, 


DC DCG absohito 


relative 


term 

2 ^ 2 \ T 


(lifforoucc 


(iitteroiico (%) 



(7s^/27p^/2), J = 

(7s?/2 7pi/2 7p3/2), J = 1 27198 26806 -392 -1.4 

(7s?/2 7p}/2 7pJ/2).-^ = 2 30775 30391 -384 -1.2 

(7s?/2 7P3/2),^ = 2 66068 65225 -843 -1.3 

(7si/27P3/2),-^ = 74527 73674 -853 -LI 
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TABLE V: Transition energies between different configurations of the Indium(^ = 49), Thallium(Z = 81) and Ekathallium(Z = 
113) atoms obtained by Dirac-Hartree-Fock method with Coulomb and Coulomb-Gaunt two-electron interactions. All values 
except relative difference are in cm~^. 



Indium atom 

configuration DC DCG absolute relative 

difference difference (%) 

(5s^/25p;/2) 

(5sv25P3/2) 2142 2113 -29 -1.4 

(56-2/254/2) 30436 30377 -59 -0.2 

(5s?/254/2) 30453 30393 -60 -0.2 

(5sJ/26sl/2) 22254 22201 -23 -0.2 
Thallium atom 

configuration DC DCG absolute relative 

difference difference (%) 

JSsfJ^&lhJT) 

(6S1/26P3/2) 7684 7597 -87 -1.1 

(6si/26rf3/2) 33035 32899 -136 -0.4 

(6S1/264/2) 33086 32949 -137 -0.4 

(6.S-J 2 7.sJ 2iili 21()1() -12.") -0.5 

Ekathalliurn atom 

configuration DC DCG absolute relative 

difference difference (%) 

(7sv27p;/2) 

(7sv27P:v2) 25447 25101 -346 -1.4 

(7sv274/2) 45547 45152 -395 -0.9 

(7s?/274/2) 45654 45258 -396 -0.9 

(7sv28sv2) 35336 34962 -374 -1.1 
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TABLE VI: Hyperfine structure constants in Indium(Z — 49), Thallium(Z = 81) and Ekathallium(Z — 113) atoms calculated 
by Dirac-Hartree-Fock method with Coulomb and Coulomb-Gaunt two-electron interactions for different configurations. All 
values except relative difference are in MHz. 



Indium atom 

configuration DC DCG absolute relative 
difference difference (%) 

(5s^/25pt/2) 1913 1900 ^13 IvT 

(5s?/25P3/2) 288 287 -1 -0.3 

(5s?/254/2) 4.41 4.40 -0.01 -0.2 

(5s?/254/2) 1-88 1.88 0.00 0.0 

(5s?/26sJ/2) 1013 1011 -2 -0.2 

Thallium atom 

configuration DC DCG absolute relative 

difference difference (%) 

(6si/26pi/2) 18918 18691 ^227 Tl" 

(6s?/26p3/2) 1403 1391 -12 -0.9 

(6S1/264/2) 20.8 20.8 0.0 0.0 

(63^/264/2) 8.72 8.70 -0.02 -0.2 

(6s2/27sJ/2) 7826 7807 -19 -0.2 
Ekathallium atom ^ 

configuration DC DCG absolute relative 

difference difference (%) 

(7si/27p;/2) 150168 147538 -2630 Ts 

(7s?/27pJ/2) 2007 1983 -24 -1.2 

(7si/27d3/2) 34.3 34.2 -0.1 -0.3 

(7si/274/2) 13.5 13.5 0.0 0.0 

(7s?/28si/2) 28580 28473 -107 -0.4 



The magnetic moment /ijv and spin / for the Eka-thallium nucleus were taken as those for Thallium. The presented results 
can be easily recalculated as only the proper values of and / are known because they just include the /I coefficient. 



TABLE VII: Core-core plus core-valence contributions (i.e. valence-valence contributions are excluded) to the total energy of 
the uranium atom from the Gaunt and retardation interactions for different choices of the cores (in a.u.). 





5/^7s" 


6d" 


5/^7s"6d^ 


core 


Gaunt 


Retardation 


Gaunt 


Retardation 




(DCG-DC) 


(DCB-DCG) 


(DCG-DC) 


(DCB-DCG) 


[He] 


28.62479082 


-3.07749038 


28.62643211 


-3.07782300 


[Ne] 


38.76297919 


-4.17703048 


38.76744083 


-4.17791688 


[Ar 2.d}°] 


41.53498862 


-4.40193727 


41.53928838 


-4.40254968 


[Kr 4d'"4/"] 


42.00441512 


-4.42817590 


42.00875116 


-4.42883266 


all 


42.02222003 


-4.43058931 


42.02637715 


-4.43120398 



TABLE VIII: Core-core plus core-valence contributions to the energy of the transition 5f^7s^6d^ 5/^7s^6d^ from the Gaunt 
and retardation interactions calculated from the data of Table IVlTl (in cm^'^). 



core 


[He] 


[Ne] [Ar3d"'] 


[Kr 4d^"4/^'*] 


all 


Gaunt 


360 


979 944 


952 


912 


(DCG-DC) 










Retardation 


-73 


-195 -134 


-144 


-135 


(DCB-DCG) 
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TABLE IX: The BI contributions to the total energy from different shells, calculated using eq. @ and estimated by eq. Illll 
(in cm~^). 



valence-valence contributions core-core contributions 







2 

ns 




ns^np^ 






ns^ 


atom 


BP 


o? {ns 


h\nsf-c% 


BI 


atom 


BI 


a^(ns|i|ns)^ • Cjv 


He (n=l) 


14 (14) 




19 




He (n=l) 


14 


19 


Ne (n=2) 


15 (21) 




18 


250 (488) 


Ne (n=l) 


2602 


3563 


Ar (n=3) 


3.4 (6.5) 




3.6 


55 (126) 


Ar (n=l) 


15880 


22045 


Kr (n=4) 


2.2 (5.1) 




2.2 


36 (90) 


Kr (n=l) 


133442 


197531 


Xe (n=5) 


1.3 (3.6) 




1.3 


23 (62) 


Xe (n=l) 


468959 


779730 


Rn (n=6) 


1.2 (4.1) 




1.2 


23 (64) 


Rn (n=l) 


2095503 


5170727 


U (n=7) 


0.094 (0.27) 




0.096 




U (n=l) 


2634137 


7290932 



core- valence contributions for uranium 



shells 


BI 


ns^ 

a^ns\^\nsf- 


ms^ 


=7, m=6) 


.73 




.79 


=7, m=5) 


1.9 




1.9 


=7, m=4) 


4.6 




4.3 


=7, m=3) 


11 




10 


=7, m=2) 


33 




28 


=7, m=l) 


151 




109 


=6, m=5) 


21 




18 


=6, m=4) 


51 




41 


=6, m=3) 


123 




95 


=6, m=2) 


360 




260 


=6, m=l) 


1658 
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(a) The values obtained with ea. instead of ea. ii'll are in brackets. 

(b) The normalizing coefficient, cjv — 0.34, is calculated as some average value for the BI contributions within the considered 

s-shells. 

(c) The additional factor, 2.7, to cn appears due to the fact only the exchange (Kb) terms of BI between four electron pairs 
from different shells, ns^ and ms^ (n 7^ m), are considered, whereas both the direct (Js) and exchange (Kb) BI terms 

between the only one electron pair contribute in the case of the same shell, ns^. Taking into account that \Kb \ ~ 2 ■ \ Jb\ (this 

expression is exact in the case n = m), one goes to the value 8/3 ^ 2.7. 



